Experimental evidence of a (p Josephson junction 
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We demonstrate experimentally the existence of Josephson junctions having a doubly degenerate 
ground state with an average Josephson phase ip = ±ip. The value of (p can be chosen by design 
in the interval < ip < n. The junctions used in our experiments are fabricated as 0-tt Josephson 
junctions of moderate normalized length with asymmetric and n regions. We show that (a) these 
ip Josephson junctions have two critical currents, corresponding to the escape of the phase ip from 
— p and +p states; (b) the phase ip can be set to a particular state by tuning an external magnetic 
field or (c) by using a proper bias current sweep sequence. The experimental observations are in 
agreement with previous theoretical predictions. 
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Josephson junctions (JJs) with a phase shift of it in 
the ground statefH attracted a lot of interest in re- 
cent years 2 7]. In particular, these JJs can be used 
as on-chip phase batteries for biasing various classical^ 
and quantum 0] circuits, allowing for removing exter- 
nal bias lines and reducing decoherence. Currently, it 
is possible to fabricate simultaneously both and it 
JJs using various technologies such as superconductor- 
ferromagnet heterostructures [Iol - [l^ | or JJs based on d- 
wave superconductors [l5 -18 1 . 

It would be remarkable to have a phase battery pro- 
viding an arbitrary phase shift ip, rather than just or 
7r. The simplest idea is to combine and 7r JJs to ob- 
tain a ip JJ. However, this is not as straightforward as it 
may seem. The balance between and it parts is compli- 
cated as shown in the pioneering work[19|, where condi- 
tions for ha ving a non-trivial M-state were derived. Long 



artificial [2JJ, |2l| and natural [22| - |24| arrays of . . .-0-ir-0-ir 



. . . JJs with short segments were analyzed in detail and 
suggested as systems, where a <p JJ could be realized. 
More recently, only one period of such an array, i.e^, one 
0-7r JJ, was analyzed in an external magnetic field\2ty. In 
these works, the Josephson phase <p(x) is considered as a 
sum of a constant (or slowly varying) phase ip and a de- 
viation \£(x)\ -C 1 from the average phase. Then, for the 
average phase ip one obtains an effective current-phase 
relation (CPR) for the supercurrentj2o| 
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where the averaged value of the critical current (/ c ) = 
(j c ) Lw. The CPR (pj) exactly corresponds to a p JJ 
at zero normalized magnetic field h = 0, if To < —1, cf. 
FigQJa). Here L = Lq+L^ is the total length of JJ, while 
Lq and L„ are the lengths of and it parts, accordingly, 



w is the width of J J. 

It is worth noting that the term u ip J J" , introduced 
in Ref. [20L refers to a J J with a degenerate ground state 
phase ip = ±p. In the particular case of Eq. ([l) at h = 
one has p = arccos( — l/2Fo). The coefficients Tq and T^ 
arc defined as 12 6] 



To 
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where Iq, l n are the lengths normalized to the Joseph- 
son length X,j({j c )) and j c ,o, j c ,-n are the critical current 
densities of and 7r parts, respectively. Here Xj is cal- 
culated using the average value of the critical current 
(j c ) = {LajcQ + L 7r j c , 7r )/(L + L^), 

The physics of p JJs with a CPR given by Eq. ([U is 
quite unusual |27|. I n particular, one should observe two 
critical currents 25|, |27[ at h = 0, corresponding to the 
escape of the phase from the left (—p) or the right (+p) 
well of the double- well Josephson energy potential 



U(ip) = (Ej) 



1 - cos(V') + r h hsm(ip) + sin 2 (-0) 



(3) 

where (Ej) = <I>o (I c ) /(2tt). The critical currents are 
different because the maximum slope (maximum super- 
current in Fig. Ufa)) on the rhs (positive bias) of the — p 
well is smaller than the maximum slope (maximum su- 
percurrent in Fig. [TJa)) on the rhs of the +p well, see 

Fig. nib). 

In this letter we present experimental evidences of a p 
JJ made of one and one 7r segment, see Fig. [2fa). 

The samples were fabricated as Nb|Al- 
Al 2 3 |Ni . 6 Cuo.4|Nb heterostructures [Til, HI- These 
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FIG. 1. (Color online) (a) effective CPR j a {tp) and (b) ef- 
fective Josephson energy U(ip) calculated numerically (thick 
lines) in comparison to those given by approximate ana- 
lytical formulas {TJ and ([3| (thin lines). For h < the 
U(tp) curves look mirror reflected with respect to the [/-axis, 
i.e., U(ip,h) — U(—ip,—h). (c) numerically calculated ef- 
fective Josephson energy U(ip) at h — tilted by an ap- 
plied bias current / > (energy supplied by the current 
source —&oI^)/(2tt)). Several important situations are shown: 
ground state 7 = (same curve as in (b)), I — I r (retrap- 
ping), I = 7 C _ (escape from — ip well), and I = I c + (escape 
from +Lp well). 



FIG. 2. (Color online) (a) Sketch (cross section) of the in- 
vestigated SIFS 0-7r JJ with a step in the F-layer thickness, 
(b) Optical image (top view, colored manually) of an inves- 
tigated sample having overlap geometry. The junction area 
(0-segment and 7r-segment) is 200 x 10 /im 2 . The "etched 
F-layer" area shows where the thickness of the F-layer was 
reduced from dp t7r to d,F,o to produce a 0-7T JJ. 
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FIG. 3. (Color online) Domain of existence of <p state. The * 
shows the position of the investigated JJ at T = 2.35 K. 



superconductor-insulator- ferromagnet-superconductor 
(SIFS) JJs have an overlap geometry as shown in 
Fig. [2jb). Each junction consists of two parts, a 
conventional 0- segment and a 7r-segment. It is well 
known 0, [D, HH, |3(| that the critical current in SFS or 
SIFS JJs strongly depends on the thickness dp of the 
F-layer and can become negative within some range of 
dp (it junction). Therefore, to produce the and the it 



segments, the F-layer has different thicknesses dp t Q and 
dp j7r , as shown in Fig. EJa). To achieve this, the F-layer 
of thickness cZf.tt, corresponding to j Ci7r = j c {dF l7T ) > 
(it JJ) was fabricated first. Then the area indicated in 
Fig. HJb) was etched down to g?f,0j corresponding to 
jc,o = ir. (dp t)) > (0 JJ). Usually one obtains asymetric 
j c values |3l|, i-e. j c .o 

We have studied 3 samples. One of them has very little 
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TABLE I. Junction parameters at T = 2.35 K. The values 
°f jc.o, jc,n and the normalized length L/\j({j c (x))) are ob- 
tained from the fits. 



parameter 


O-part 7r-part 


whole JJ 


physical length 


100 ^m 100 ^m 


200 [im 


jc,o, jc,K, (jc(x)} (A/cm 2 ) 


+62.9 -47.9 


+7.5 


normalized lengths 


0.68 0.68 


1.36 



j c asymmetry of about 5%, and Lq = L n = 25 /xm, which 
corresponds to Lo/Aj(jc,o) ~ L v / 'Xj(jc^) ~ 0.37. In this 
case even 5% of asymmetry brings the sample out of the 
•p domain, see Fig. [3] for a qualitative picture. The other 
two JJs have Lq = = 100 /im and are deep inside the 
(^-domain in parameter space. Both samples show similar 
results. Here we present the results obtained on one of 
them. Its parameters are summarized in Tab. U and its 
position within the tp domain is indicated in Fig. [31 The 
values of j Cl o and j Ci „ cannot be measured directly. In 
our case we have measured the I C {H) dependence (see 
below) and then simulated it numerically using and 
j Cj7r as fitting parameters. The best fitting was obtained 
for the values specified in Tab. |TJ 

According to theoretical predictions 2^, 27] for a tp J J 
at zero magnetic field H = one expects two critical 
currents |J C _| < |J C +| (for each bias direction), corre- 
sponding to the escape of the phase from —ip and +tp 
wells of U(ip) , respectively for bias current 7 > and 
vice versa for / < 0. The current I c + is always observed. 
To observe I c - one has to have low damping so that re- 
trapping in the +tp well is avoided (for positive I). The 
I—V characteristic (IVC) of the investigated 0-ir JJ at 
H = measured at T ~ 4.2 K show only one critical 
current. Therefore the experiments were performed in a 
3 He cryostat at T down to 300 mK where the damping 
in SIFS JJs reduces drasticallvflil]. In the temperature 
range from 3.5 K down to 300 mK both I c+ and I c - are 
clearly visible in the IVCs as shown in Fig. 21 

Earlier [2?} we proposed a technique that allows to 
choose which critical current one traces in the IVC, i.e. 
from which well the phase escapes. The control is done 
by choosing a proper bias sweep sequence. For exam- 
ple, if the junction is returning from the positive voltage 
state, the potential U(tfj) is tilted so that the phase slides 
to the right. When the tilt becomes small enough the 
phase will be trapped, presumably in the right +tp well, 
cf. Fig. [He). However, this natural assumption is not 
always true (see below) . Then, if the phase is trapped at 
+tp, we can sweep the bias (a) in the positive direction 
and will observe escape from +ip (to the right) at I c+ or 
(b) in the negative direction to observe escape from +^p 
(to the left) at —I c -- 

In experiment, at T < 2.3 K, when the damping is 
very low, the currents ±7 C + and ±7 C _ are traced in ran- 
dom order. Recording one IVC after the other, we were 
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FIG. 4. (Color online) Current-voltage characteristics show- 
ing lower ±7 C _ and higher ±J C + critical currents measured 
at T m 2.35 K. At this temperature the behavior is deter- 
ministic: if one sweeps / as — J max — > +Imax -/max one 
always observes the critical currents ±/ c _; if one sweeps / as 
+/max — > — > +/ m ax, one always observes +/ c +; finally, if 
one sweeps I as — J max — ► — > —/max, one always observes 

-Ic+- 



able to obtain IVCs with all 4 possible combinations: 
(a) (-J c _,+J c _), (b) (-I c _,+I c+ ), (c) (-J c+ ,+/ c _), (d) 
(— 7 C +, +/ c +). Choosing a specific sweep sequence as de- 
scribed above does not make the outcome (I c — or I c +) 
predictable. We believe that in this temperature range 
the damping is so low that, upon returning from the pos- 
itive voltage state, the phase does not simply stop in the 
+ip well, but can also reflect from the barrier and find 
itself in a — cp well, cf. Fig. []Jc). The absence of deter- 
minism suggests that most probably we are dealing with 
a system exhibiting chaotic dynamics. This issue will be 
investigated elsewhere. Nevertheless, at T w 2.35 K we 
managed to achieve deterministic behavior as described 
above, see Fig. 2] 

Another fingerprint of a ip J J is its I C (H) dependence, 
which (a) should have the main cusp-like minima shifted 
off from H = point-symmetrically with respect to the 
origin, see Fig. 4 of Ref. [2^; (b) should show up to four 
branches in total (for both sweep polarities) at low mag- 
netic field 25j. In essence, the latter feature results from 
the escape of the phase from two different energy min- 
ima in two different directions and is an extension of the 
two-critical-currents story to the case of non-zero mag- 
netic field. Instead, the feature (a) alone cannot serve 
as a proof of a ip JJ as it is a common feature of every 
asymmetric 0-tt J J even if its ground state is or 7r[2oT] . 

The experimentally obtained ±I C (H) dependence at 
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FIG. 5. (Color online) Experimentally measured depen- 
dence Ic(H) (black symbols) and numerically calculated curve 
(red/gray smaller symbols) that provides the best fit. 



T = 2.35 K is shown in Fig. [SJ First, in the whole 
temperature range 0.3-4.2 K we observe the main min- 
ima shifted point-symmetrically to a finite field, which is 
H min = Tj ■ Icoii = <q • ±879 fiA. Second, at T < 3.5 K 
one observes the crossing of two branches and two criti- 
cal currents for each bias current polarity. The left hand 
branches correspond to the escape of the phase out of the 
+</?-well, while the right branch corresponds to the escape 
from the — (p-we\\. Note, that at T = 4.2 K the crossing 
of branches is not visible. One just observes a cusp-like 
minimum. In this case, both I c ± cannot be seen together 
and the existence of two states, as should be present for 
a ip-33, cannot be proven. 

To trace the intersection of the branches better we were 
applying a special value H reset of the magnetic field dur- 
ing reset from the voltage V > state back to V = 
state, to "prepare" a specific state (+cp or — (p) of the sys- 
tem. Then the field was set back to the "current" value 
H and the bias was ramped up to trace I C (H). By doing 
this for different values of H leset , one is able to trace the 
intersection of the branches much better than without 
this technique, see Fig. [5l Still, in experiment we were 
not able to trace the branches for positive and negative 
I c up to the point where they meet, as shown in Fig. 4(a) 
and (b) of Ref. [25J, probably because of retrapping. 

To extract some parameters of our J J from the I C {H) 
curve, we performed numerical simulations, by solving 
the full sine-Gordon equation for a 0-7T JJ, using the nor- 
malized length I and the critical currents j Ci o and j Cj7r as 
fitting parameters. The simulations were performed us- 
ing STKjj[32j]. The objective was to obtain the best 
fit close to the origin, especially the point where the 
branches of I C {H) cross. The best fit is shown in Fig. [5] 
and was obtained for I — 3.7 and the values jc,o and jc,<K 
from Tab. [TJ Measured data in Fig. [5] are shifted along 



I co ii axis by 530 /iA to compensate for average remanent 
magnetization of the F-layer. To obtain an almost per- 
fect fit we have also assumed a difference in the constant 
remanent magnetizations in the and in the 7r parts of 
A (M) = ?7l.5mA, similar to earlier studies (28. 31 1. One 
can also see that one of the experimental branches after 
the main maximum runs parallel to the simulated curve 
(the experimental I C (H) is not point symmetric). This 
shift stays the same even if we cycle the sample through 
T c of Nb, but changes if we cycle through the Curie tem- 
perature Tc of the F-layer. We conclude that this is re- 
lated to the non-uniform remanent magnetization of the 
F-layer. 

The parameters obtained from the fit allow to see the 
location of our JJ on the (lo,lir) plane, see Fig. [3l One 
sees that the JJ is not really short and lays quite deep in 
the if domain. In this region the analytic results [25, 26 1 
and, in particular, Eqs. ([1]) and ([3]) are valid only qual- 
itatively, c/., Fig. 3 in Ref. 0. Therefore we calculated 
the CPR of our JJ numerically. We started with a 0-7T JJ 
with parameters Iq, l w , j Cj o and j„ obtained from the fit, 
assumed some applied bias current / and found all static 
solutions (j)(x) numerically. Then for each of these solu- 
tions we calculated tp = (4>(x)) and plotted all those i/j 
on a plot. By repeating this procedure for different 
/, we obtain an effective CPR shown in Fig.QJa) to- 
gether with the curve produced by the analytical formula 
(HJ). One can see that the exact effective CPR calculated 
numerically and the approximate CPR calculated ana- 
lytically are qualitatively similar. From the numerical 
CPR the value of <p = OAln for our particular (p JJ is 
extracted. 

In conclusion, we have demonstrated experimentally 
the realization of a ip Josephson junction based on a 0-7r 
SIFS junction. We have observed experimentally a shift 
of the I C {H) minimum according to the effective CPR 
([T]), as predicted [2^], as well as two critical currents ±/ c ± 
close to zero magnetic field for each bias current polarity. 
We showed that one can choose between the ±tp states 
of the system using an externally applied magnetic field, 
which removes their degeneracy. We also showed that one 
can bring the system into one of the two states ±ip by 
properly tilting the potential using the bias current. De- 
pending on the damping (temperature) one can achieve 
deterministic behavior (when damping is small enough 
to see the lower critical current, but large enough to trap 
the phase in a particular well) as well as random behav- 
ior at very low damping. The obtained JJ has tp = 0.417T 
at 1= 0). 
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